Retarded coordinates based at a world line, and the motion of a small black hole in an 

external universe 
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In the first part of this article I present a system of retarded coordinates based at an arbitrary 
world line of an arbitrary curved spacetime. The retarded-time coordinate labels forward light 
cones that are centered on the world line, the radial coordinate is an affine parameter on the null 
generators of these light cones, and the angular coordinates are constant on each of these generators. 
The spacetime metric in the retarded coordinates is displayed as an expansion in powers of the 
radial coordinate and expressed in terms of the world line's acceleration vector and the spacetime's 
Riemann tensor evaluated at the world line. The formalism is illustrated in two examples, the first 
involving a comoving world line of a spatially-flat cosmology, the other featuring an observer in 
circular motion in the Schwarzschild spacetime. The main application of the formalism is presented 
in the second part of the article, in which I consider the motion of a small black hole in an empty 
external universe. I use the retarded coordinates to construct the metric of the small black hole 
perturbed by the tidal field of the external universe, and the metric of the external universe perturbed 
by the presence of the black hole. Matching these metrics produces the MiSaTaQuWa equations of 
motion for the small black hole. 



I. INTRODUCTION 

In the first part of this article I present a system of re- 
tarded coordinates (u, r, 8, (f>) based at an arbitrary world 
line 7 of an arbitrary spacetime with metric g a p. The co- 
ordinates are adapted to the forward light cone of each 
point z(t) of the world line (r is the proper-time param- 
eter on 7). The retarded-time coordinate u is constant 
on each light cone, and it agrees with r at the cone's 
apex. The radial coordinate r is an affine parameter on 
the cone's null generators, and it gives a measure of dis- 
tance away from the world line. The angular coordinates 
6 A = (9, <f>) are constant on each of these generators. 

This formalism complements a line of research that 
was initiated by Synge [ll and pursued by Ellis and his 
collaborators Q, 01 HI TH m their work on observa- 
tional cosmology. While the central ideas exploited here 
are the same as with Synge and Ellis, my implementation 
is substantially different: While Synge and Ellis sought 
definitions for their optical or observational coordinates 
that apply to large regions of the spacetime, my consid- 
erations are limited to a small neighborhood of the world 
line; and while Ellis favored a construction based on past 
light cones (because of the cosmological context in which 
information at 7 is gathered from the past), here the 
preference is given to future light cones. 

The introduction of retarded coordinates is motivated 
by the desire to construct solutions of wave equations 
for massless fields that are produced by pointlike sources 
moving on the world line. The retarded coordinates natu- 
rally incorporate the causal relation that exists between 
the source and the field, and for this reason the solu- 
tion takes a simple explicit form (in the neighborhood in 
which the coordinates are defined). The point of view 
developed here, therefore, is opposite to the cosmological 
view developed by Ellis: here 7 is the cause and the effect 



is propagated along forward light cones. 

A quasi-Cartesian version of the retarded coordinates 
is introduced in Sec. II B, after reviewing some necessary 
geometrical elements (a tetrad transported on 7, Synge's 
world function, and the bitensor of parallel transport) in 
Sec. II A. In Sec. II C, I explore the significance of r as 
an affine parameter on the generators of the light cones, 
and describe the vector k a that is tangent to the congru- 
ence of generators. The metric in retarded coordinates 
is gradually constructed in Sees. II D through F, and its 
quasi-Cartesian form is displayed in Eqs. I|2.39[l - i|2.41|l . 
In Sec. II G, I carry out a transformation to angular 
coordinates, and this form of the metric is displayed in 
Eqs. <|2.55f) - H2.58fl . I n Sec. II H, I present an impor- 
tant simplification of the formalism that occurs when 7 
is a geodesic and the Ricci tensor vanishes on the world 
line. Two examples are worked out in Sees. II I and J: 
first I apply the formalism to a comoving world line of a 
spatially-flat cosmology, and then I examine the metric 
near an observer in circular motion in the Schwarzschild 
spacetime. 

The main application of the formalism is presented in 
the second part of the paper. Here I consider the mo- 
tion of a nonrotating black hole of small mass m in a 
background spacetime with metric g a p; the metric is as- 
sumed to be a solution to the Einstein field equations in 
vacuum. By employing the powerful method of matched 
asymptotic expansions (see, for example, Ref. |8| for a 
discussion), I derive equations of motion for the black 
hole. If z m (t) are the parametric relations describing the 



black hole's world line, 



dz^/dr is the velocity vec- 



tor, and a M = Du^/dr the acceleration vector, then the 
equations of motion take the form of 
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where 

x/'/dr' (1.2) 

is an integration over the past portion of the world 
line; the integral involves the retarded Green's function 
G^(j.'v'(z(t), z(t')) for the metric perturbation |j| as- 
sociated with the small black hole, and it is cut off at 
t' = t~ = t — + to avoid the singular behavior of the 
Green's function at coincidence. 

Equation l|l.l|l is not new, and the method of deriva- 
tion presented here also is not new; they both originate 
in a 1997 paper by Mino, Sasaki, and Tanaka [T^j. (The 
equations were later rederived by Quinn and Wald |ll| . 
and they are now known as the MiSaTaQuWa equa- 
tions of motion.) The method of matched asymptotic 
expansions was first used to establish that in the limit 
in — > 0, the motion of a black hole is geodesic in the 
background spacetime 0, 0, [w| . Using this approach, 
Mino, Sasaki, and Tanaka (see also Ref. ^5|) were able 
to compute the first-order correction to this result, and 
showed that at order m, the motion is accelerated and 
governed by Eq. i|l.lfl . Their derivation, however, did 
not rely on a specific choice of coordinates, and it left 
many details obscured. By adopting retarded coordi- 
nates, these details can be revealed and the matching 
procedure clarified. This was my aim here: to provide 
a clearer, more explicit implementation of the matching 
procedure employed by Mino, Sasaki, and Tanaka. 

The method of matched asymptotic expansions is ex- 
plained in Sec. Ill A. Essentially, the metric of the small 
black hole perturbed by the tidal gravitational field of the 
external universe is matched to the metric of the external 
universe perturbed by the black hole; ensuring that this 
metric is a valid solution to the vacuum Einstein field 
equations determines the motion of the black hole. The 
perturbed metric of the small black hole is constructed 
in an internal zone in terms of internal retarded coor- 
dinates (u,i,6 A ); this is carried out in Sec. Ill B. The 
perturbed metric of the external universe is constructed 
in an external zone in terms of external retarded coor- 
dinates (u,r,9 A ); this is carried out in Sec. Ill C. The 
matching of the two metrics is performed in a buffer zone 
that overlaps with both the internal and external zones; 
this involves a transformation from the external to the 
internal coordinates that is explicitly worked out in Ap- 
pendix C. The matching is carried out in Sec. Ill D, where 
we see that it does indeed produce Eq. I|l.lfl . 

This concludes the summary of the results contained 
in this paper. Various technical details are relegated to 
Appendices A and B. Throughout the paper I work in 
geometrized units (G = c = 1) and with the conventions 
of Misner, Thorne, and Wheeler ^(|. Other applications 
of the retarded coordinates (including the fields produced 
by point particles, and their motion in curved spacetime) 
are featured in a recent review article 1171 . 



II. RETARDED COORDINATES 
A. Geometrical elements 

To construct the retarded coordinates we must first in- 
troduce some geometrical elements on the world line 7 at 
which the coordinates are based. The world line is de- 
scribed by relations z^(t) with r denoting proper time, 
its normalized tangent vector is = dz^/dr, and its 
acceleration vector is = Du^/dr with D/dr denoting 
covariant differentiation along the world line. Through- 
out we use Greek indices fi, v 1 A, p, etc. to refer to tensor 
fields defined, or evaluated, on the world line. 

We install on 7 an orthonormal tetrad that consists of 
the tangent vector u M and three spatial vectors e^. These 
are transported on the world line according to 

i^L=a a u» + u> a b el (2.1) 

where a a (r) — a^e^ are the frame components of the ac- 
celeration vector and u! a b(T) = — Wfeo(r) is a prescribed 
rotation tensor. Setting uj a b — would make the triad 
Fermi- Walker transported on the world line, and in many 
applications this would be a sensible choice. We never- 
theless allow the spatial vectors to rotate as they are 
transported on the world line; this makes the formalism 
more general, and we will need this flexibility in Sec. Ill 
of the paper. It is easy to check that Eq. I|2.1|l is com- 
patible with the requirement that the tetrad (w^e^) be 
orthonormal everywhere on 7. 

From the tetrad on 7 we define a dual tetrad (e°,e^) 
with the relations 

e° - e a ^ = 5 ah 9tlv el. (2.2) 

The dual vectors eJJ satisfy a transport law that is very 
similar to Eq. J23J- The tetrad and its dual give rise to 
the completeness relations 

(2.3) 

9nv = ~ e fj, e ii $ab 

for the metric and its inverse evaluated on the world line. 

The retarded coordinates are constructed with the help 
of a null geodesic that links a given point x to the world 
line. This geodesic must be unique, and we thus restrict 
x to be within the normal convex neighborhood of 7. We 
denote by (3 the unique, future-directed null geodesic that 
goes from the world line to x, and x' = z{u) is the point 
at which f) intersects the world line; u is the value of the 
proper-time parameter at this point. To tensors at x we 
shall assign the Greek indices a, (3, 7, 6, etc.; to tensors 
at x' we shall assign the indices a' , /?', 7', 6', and so on. 

From the tetrad (u a , e" ) at x 1 we construct another 
tetrad (eQ,e") at x by parallel transport on (3. By rais- 
ing the frame index and lowering the vectorial index, 
we obtain also a dual tetrad at x: = — gape® and 
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e a = ^ ab ffa/3 e b • The metric at x can then be expressed 
as 

g a p = -e°e° + 5 ah e a a e% (2.4) 

and the parallel propagator yjj (also known as the bivec- 
tor of geodetic parallel displacement 18]) from x' to x is 
given by 

g a a ,{x,x r ) = -e%u a > + e%e a a , , 

(2.5) 

g a ' a {x\x) = u a 'e° a + e:'e a a . 

This is defined such that if A a is a vector that is parallel 
transported on (3, then A a (x) = g a a ,(x,x')A a (x') and 
A a (x') — g a a {x\ x)A a {x). Similarly, if p a is a parallel- 
transported dual vector, then p a (x) = g a a (x' , x)p a i{x') 
and p a >{x') = g a a ,(x,x')p a (x). 

The last ingredient we shall need is Synge's world func- 
tion a (z, x ) jl| (also known as the biscalar of geodetic in- 
terval [l8|)- This is defined as half the squared geodesic 
distance between the world-line point z(r) and a neigh- 
boring point x. The derivative of the world function with 
respect to z^ is denoted <Jp{z, x); this is a vector at z (and 
a scalar at x) that is known to be tangent to the geodesic 
linking z and x. The derivative of cr(z, x) with respect to 
x a is denoted a a (z,x); this vector at x (and scalar at z) 
is also tangent to the geodesic. We use a similar notation 
for multiple derivatives; for example, a^ a = Vq-V^ct and 
<7 a fj = V^VqCt, where V a denotes a covariant derivative 
at x while V M indicates covariant differentiation at z. 

The vector —a^(z, x) can be thought of as a separation 
vector between x and z, pointing from the world line to 
x. When x is close to 7, —a^(z,x) is small and can be 
used to express bitensors in terms of ordinary tensors at 
z 0,0. For example, 

<?nu = g^u - ^R^\u P <y x <J p H — , (2.6) 
^ a = -0q + ^R^u P <J x <J p + ■ ■ •) , (2.7) 

v a p = 9^9 U (g^ - \r^ p ^<t p + •••), (2.8) 

where g^ a = g^ a (z,x) is the parallel propagator and 
Rp\vp is the Riemann tensor evaluated on the world line. 



B. Definition of the retarded coordinates 



The last statement indicates that x' = z(u) and x are 
linked by a null geodesic. 1 

From the fact that a a is a null vector we obtain 

r=(S ab x a x b )^ 2 =u a >a a ', (2.10) 

and r is a positive quantity by virtue of the fact that (3 
is a future-directed null geodesic — this makes a a past- 
directed. In flat spacetime, a a = —(x — x') a , and in 
a Lorentz frame that is momentarily comoving with the 
world line, r = t — t' > 0; with the speed of light set equal 
to unity, r is also the spatial distance between x' and x as 
measured in this frame. This gives us an interpretation 
of r = u a :<j a as a retarded distance between x and the 
world line, and we shall keep this interpretation even in 
curved spacetime. The claim that r gives a measure of 
distance between x' and x will be substantiated in Sec. II 
C. 

Another consequence of Eq. (|2.9fl is that 

a a ' = -r(u a ' +n a e*'), (2.11) 
where ft a = x a /r is a frame vector that satisfies 

s ab n a n b = 1. 

A straightforward calculation reveals that under a dis- 
placement of the point x, the retarded coordinates change 
according to 

du = -k a dx a , (2.12) 
dx a = ~(ra a ~uj a b x b + e a a ,cryu l3 ')du 

-e a a ,a a ' p dx^, (2.13) 

where k a = a a jr is a future-directed null vector at x 
that is tangent to the geodesic j3. To obtain these results 
we must keep in mind that a displacement of x typically 
induces a simultaneous displacement of x', because the 
new points x + Sx and x' + Sx' must also be linked by a 
null geodesic. We therefore have = a(x + 5x, x' + Sx') = 
a a 5x a + cr a > 8x a , and Eq. I|2.12|l follows from the fact 
that a displacement along the world line is described by 
8x a = u a Su. 



C. Retarded distance; null vector field 

If we keep x' linked to x by the relation a{x, x 1 ) = 0, 
then 

r{x) =a a ,(x,x')u a '(x') (2.14) 

can be viewed as an ordinary scalar field defined in a 
neighborhood of 7. We can compute the gradient of r by 



The quasi-Cartesian version of the retarded coordi- 
nates are defined by 

x° = u, x a — — e a a ,(x)a a (x,x), a(x,x') = 0. 

(2.9) 



A s imilar definition can be given for Fermi normal coordinates Qj, 
lUl . Here x is linked to 7 by a spacelike geodesic that intersects 
the world line orthogonally. The intersection point is such that 
<Jp(x, z)u^ = 0, and this replaces the last condition of Eq. 12.91 . 
The other relations are unchanged. 
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finding how r changes under a displacement of x (which 
again induces a displacement of x'). The result is 

dpr — -(ova Q + a a 'fi'U a vP )kp + a a 'pu a . (2.15) 
Similarly, we can view 

*■(*> = 

as an ordinary vector field, which is tangent to the con- 
gruence of null geodesies that emanate from x' . It is easy 
to check that this vector satisfies the identities 

a a phP = kc, a a > kP = (2.17) 

from which we also obtain a a 'pu a k^ = 1. From this last 
result and Eq. (|2.15|) we deduce the important relation 

k a d a r = 1. (2.18) 

In addition, combining the general statement a a — 
-g a a ,(T a ' with Eq. (|2TT|) gives 

k a =g a a ,(u a ' +n a e%'); (2.19) 

the vector at x is therefore obtained by parallel transport 
of u a + O a e" on 0. From this and Eq. (|2.4|l we get the 
alternative expression 

k a = e« + fl a e«, (2.20) 

which confirms that k a is a future-directed null vector 
field (recall that f2 a = x a jr is a unit frame vector). 

The covariant derivative of k a can be computed by 
finding how the vector changes under a displacement of 
x. (It is in fact easier to first calculate how rk a changes, 
and then substitute our previous expression for dpr.) The 
result is 

rk a ;[j = U a p — fc a tT/3 7 'U 7 — kfiOa^'U 1 

+ {a a >a a ' + a a 'p>u a 'u ')k a kp. (2.21) 

From this we infer that k a satisfies the geodesic equa- 
tion in affine-parameter form, k^.ak^ = 0, and Eq. (|2.18|1 
informs us that the affine parameter is in fact r. A dis- 
placement along a member of the congruence is therefore 
described by dx a — k a dr. Specializing to retarded coor- 
dinates, and using Eqs. I|2.12[) . (|2.13l) . and (|2.17l) . we find 
that this statement becomes du = and dx a = (x a /r) dr, 
which integrate to u — constant and x a — rfl a , respec- 
tively, with n a representing a constant unit vector. We 
have found that the congruence of null geodesies emanat- 
ing from x' is described by 

u = constant, x a = rtt a {0 A ) (2.22) 

in the retarded coordinates. Here, the two angles 9 A 
(A = 1,2) serve to parameterize the unit vector fl a , 
which is independent of r. 



Equation l|2.21|l also implies that the expansion of the 
congruence is given by 

n a — 2 

6 = fc" = — ^ . (2.23) 

r 

Using the expansion for a a a given by Eq. (|2.8|) . we find 
that this becomes r9 = 2 - \R a <po a 'a 13 ' + 0(r 3 ), or 

rO = 2- ir 2 (Rao + 2R aa Q a + R ab n a n b ) +0(r 3 ) (2.24) 

after using Eq. I|2.11|l . Here, i?oo = R a 'P'U a u 13 , Ro a = 

R a '(3'U a , and R a b — Ra'p'e" are the frame com- 
ponents of the Ricci tensor evaluated at x'. This result 
confirms that the congruence is singular at r = 0, because 
diverges as 2/r in this limit; the caustic coincides with 
the point x' . 

Finally, we infer from Eq. (|2.21|) that k a is hypersur- 
face orthogonal. This, together with the property that k a 
satisfies the geodesic equation in affine-parameter form, 
implies that there exists a scalar field u(x) such that 

k a = -d a u. (2.25) 

This scalar field was already identified in Eq. I|2.12|l : it 
is numerically equal to the proper-time parameter of the 
world line at x' . We may thus conclude that the geodesies 
to which k a is tangent are the generators of the null cone 
u = constant. As Eq. I|2.22|l indicates, a specific genera- 
tor is selected by choosing a direction 11° (which can be 
parameterized by two angles 9 A ), and r is an affine pa- 
rameter on each generator. The geometrical meaning of 
the retarded coordinates is now completely clear, and r is 
recognized as a meaningful measure of distance between 
x and the world line. The construction is illustrated in 
Fig. 1. 

D. Frame components of tensors on the world line 

The metric at x in the retarded coordinates will be 
expressed in terms of frame components of vectors and 
tensors evaluated on 7. For example, if a a is the accel- 
eration vector at x' , then as we have seen, 

«a(«)=ve:' (2.26) 

are the frame components of the acceleration at proper 
time u. 
Similarly, 



RaObo(u) 


= Ra 


'Yf3>S> 


e Q ' 


U 1 ' 


' (3' 5' 


RaObd(u) 


= Ra 


'Y/3'S' 


e a ' 


U 1 ' 


e b e d 


Racbd(u) 


= Ra 




e a ' 


e 1 ' 


f3' 6' 
e b e d 



are the frame components of the Riemann tensor evalu- 
ated on 7. From these we form the useful combinations 

Sab = RaObO + RaObc^ C + RbOac^ C + Racbd^ C ^ d 
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FIG. 1: Retarded coordinates of a point x relative to a world 
line 7. The retarded time u selects a particular null cone, 
the unit vector Q a = x a /r selects a particular generator of 
this null cone, and the retarded distance r selects a particular 
point on this generator. 



These results can also be expressed in the form of gra- 
dients of the retarded coordinates: 



d a u = e° - £l a e a a , 



(2.34) 



d a x a = 



1 



ruffi + -r^S a + 0{r 3 ) 



(ra a 



2 

roj a r tt c 



r 2 s a )n b 



+ -r"S a b + 0(r 3 ) 



(2.35) 



Notice that Eq. (|2.34|) follows immediately from 
Eqs. l{2~2"0"|) and (|2~2"5|) . From Eq. l(235|l and the iden- 
tity d a r = fl a d a x a we also infer 

ra a n a + l-r 2 S + O(r 3 )}e a 



d a r = - 



1 



2 

ra b tt b 



-r 2 S a + 0(r 3 ) 




(2.36) 



where we have used the facts that S a = S ab fl b and S — 
S a n a ; see Eqs. QFZfy and l(2~3Uj). It may be checked that 
Eq. l(2~3tJj) agrees with Eq. lt2~T5l) . 



= S ba , (2.28) 

S a = S ab Q b = RaObO^ ~ RabOc^^, (2.29) 

S = S a Q a = R a o bQ n a Q b , (2.30) 

in which the quantities tt a = x a /r depend on the angles 
9 A only — they are independent of u and r. 

We have previously introduced the frame components 
of the Ricci tensor in Eq. (|2.24|) . The identity 

R 00 + 2R Qa n a + R ab n a n b = s ab s ab - s (2.31) 



follows easily from Eqs. 
of the Ricci tensor. 



and the definition 



E. Coordinate displacements near 7 

The changes in the quasi-Cartesian retarded coordi- 
nates under a displacement of x are given by Eqs. (|2.12|) 
and (|2.13|l . In these we substitute the expansions for 
<j a 'f)' an d a a' (3 that appear in Eqs. (|2.6|l and 12. 7f) . as 
well as Eqs. (|2.11|l and l|2.19|l . After a straightforward 
calculation, we obtain the following expressions for the 
coordinate displacements: 

{e° a dx a )-Q a (e b a dx a ), 
1 



du 
•is" = 



ra ~ rto b il 



„2 S a 



0{r 3 ) 

5 a b + (ra a - ruj a c il c + K 2 S a ^Q 
+ \r 2 S\ + <3(r 3 )] (e b a dx a ). 



(2.32) 

s° dx a ) 



(2.33) 



Notice that the result for du is exact, but that dx a is 
expressed as an expansion in powers of r. 



F. Metric near 7 

It is straightforward to invert the relations of 
Eqs. (|2~32l and iflHSj) and solve for e° dx a and e a a dx a . 
The results are 



l + ra a n a + ^r 2 S + 0{r 3 ) 



du 



l + -r 2 S)<J„ 



(> 



r 2 S a +0{r 3 ) 



dx a , 



r(a a - uj a b n b ) + -r 2 S a + 0{r 3 ) 



+ -r z S a n b + 0(r 3 ) 



dx" 



(2.37) 
du 

(2.38) 



These relations, when specialized to the retarded coordi- 
nates, give us the components of the dual tetrad (e° , e°) 
at x. The metric is then computed by using the com- 
pleteness relations of Eq. (|2.3|l . We find 



- (1 + ra a il a ) 2 + r 2 (a a - uj ab n b ) (a a ~ <ft c ) 
-r 2 S + 0(r 3 ), (2.39) 



9ab 



= -yl + ra b il + -r 2 SJQ a + r(a a - iv ab Q 
+ r 3 r 2 S a + 0(r 3 ), 

i + ^r 2 s)n a n b -^r 2 s ab 

+ ^r 2 (S a n b + n a S b ) +0(r 3 ). 



(2.40) 



(2.41) 
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We see that the metric possesses a directional ambiguity 
on the world line: the metric at r = still depends on the 
vector il a — x a jr that specifies the direction to the point 
x. The retarded coordinates are therefore singular on the 
world line, and tensor components cannot be defined on 
7. Because we are working with frame components of 
tensors, this poses no particular difficulty. 

By setting S ab = S a = S = in Eqs. ljZ55 |l -l|^lT )l 
we obtain the metric of flat spacetime in the retarded 
coordinates. This we express as 

Vuu = -(1 + ra a n a ) 2 +r 2 (a a -iu ab n b )(a a -Lu a c Q c ), 
Vua = -(l + ra b n b )n a + r(a a -uj ab n b ), (2.42) 

Vab = <>ab — ^a^b, 

and we see that the directional ambiguity persists. This 
should not come as a surprise: the ambiguity is present 
even when a a — w a b — 0, and is generated simply 
by performing the coordinate transformation u = t — 
\/ x 1 + y 2 + z 2 . The retarded coordinates are therefore 
necessarily singular at the world line. But in spite of the 
directional ambiguity, the metric of flat spacetime has a 
unit determinant everywhere, and it is easily inverted: 

f] uu = 0, 

7 f a = -n a , (2.43) 
n ab = 5 ab + r(a a -u> a c n c )n b + rfl a (a b ~u> b c n c ). 

The inverse metric also is ambiguous on the world line. 

To invert the curved-spacetime metric of Eqs. I|2.39[) - 
(|2.41|) we express it as g a p = i] a p + h a p + 0(r 3 ) and 
treat h a p = 0(r 2 ) as a perturbation. The inverse metric 
is then g afi = if? - n^rj^h^s + 0(r 3 ), or 



0. 



,. ab 



-n a , 



(2.44) 
(2.45) 



= S ab + r(a a -Lu a c Q c )n b + rSl a (a l 
+ \r 2 S ab + \r 2 (S a n b + Q a S b ) + 0(r 3 ). (2.46) 

The results for g uu and g ua are exact, and they fol- 
low from the general relations g a ^{d a u){dpu) = and 
g a,3 (d a u){dpr) = —1 that are derived from Eqs. (|2.18|) 
and (Egg) . 

The metric determinant is computed from \J—g = 1 + 
\rf h a(i + 0(r 3 ), which gives 



l--r 2 (5 ab S ab -S)+0(r 3 ) 

1 - ir 2 (i? 00 + 2R 0a n a + Ra b n a n b ^ 

+ 0(r 3 ), 



(2.47) 



where we have substituted the identity of Eq. I|2.31|l . 
Comparison with Eq. (|2.24|l gives us the interesting re- 
lation y/— g = ir0 + 0(r 3 ), where is the expansion of 
the generators of the null cones u = constant. 



G. Transformation to angular coordinates 



Because the frame vector fl a = x a /r satisfies 
5 ab il a £l b — 1, it can be parameterized by two angles 
9 A . A canonical choice for the parameterization is fl a = 
(sin cos 4>, sin 9 sin cf>, cos 9) . It is then convenient to per- 
form a coordinate transformation from x a to (r, 9 A ), us- 
ing the relations x a — rVl a {9 ). (Recall from Sec. II C 
that the angles 9 A are constant on the generators of the 
null cones u = constant, and that r is an affine parame- 
ter on these generators. The relations x a = rf2 a therefore 
describe the behavior of the generators.) The differential 
form of the coordinate transformation is 



dx a = fl a dr + rn a A d9 A , 
where the transformation matrix 

ilA ^d9^ 
satisfies the identity Vl a U, a A = 0. 
We introduce the quantities 

&ab = 5 ab Vt a A Vt b B , 



(2.48) 



(2.49) 



(2.50) 

which act as a (nonphysical) metric in the subspace 
spanned by the angular coordinates. In the canonical 
parameterization, Qab = diag(l, sin 2 9). We use the in- 
verse of Qab, denoted Sl AB , to raise upper-case Latin 
indices. We then define the new object 

(2.51) 



5 a bQ AB Q b R 



which satisfies the identities 



n A n% 



s A 



o a o A 



5\ - n a n b . 



(2.52) 



The second result follows from the fact that both sides 
are symmetric in a and b, orthogonal to fl a and il b , and 
have the same trace. 

From the preceding results we establish that the trans- 
formation from x a to (r, 6 ) is accomplished by 



dx a 
dr 



d9 A ~ rilA > 



(2.53) 



while the transformation from (r,9 A ) to x a is accom- 
plished by 



— = n«, 

dx a 



— = v. 

dx a r a 



(2.54) 



With these rules it is easy to show that in the angular 
coordinates, the metric takes the form of 

ds 2 = g uu du 2 ~ 2 dudr + 2g uA dud9 A + g AB d9 A d9 B , 

(2.55) 

with 



9uu 

9uA = 
9AB = 



-(1 + ra a n a ) 2 + r 2 (a a - uj ab n b ) (a a - u a ,Sl c ) 
-r 2 S + 0(r 3 ), (2.56) 



r(a a - to ab n b ) + -r 2 S a + 0(r 3 )\ Q% (2.57) 



n AB - -r 2 s ab n a A n b B + 0(r 3 ) 



(2.58) 
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The results g ru = — 1, g rr = 0, and g r A — are ex- 
act, and they follow from the fact that in the retarded 
coordinates, k a dx a = —du and k a d a = d r . 

The nonvanishing components of the inverse metric are 



9 ur = -1, 

g rr = l + 2ra a n a + r 2 S + 0(r 3 ), 



,rA 



9 



AB 



1 - 

r - 
+ 0{r 3 ) 
1 



~(a a -io a b n b ) + ^r 2 S a 



(2.59) 
(2.60) 



(2.61) 
(2.62) 



The results g uu = 0, g ur = — 1, and g uA — are exact, 
and they follow from the same reasoning as before. 

Finally, we note that in the angular coordinates, the 
metric determinant is given by 



-g = r 



l-\r 2 {Roo+2Roatt a +R ab n a n b ) +0(r 3 ) 



(2.63) 

where f2 is the determinant of Qab', in the canonical 
parameterization, — smc/. 



H. Specialization to o* 1 = = Rp, v 

In this subsection we specialize our previous results to 
a situation where 7 is a geodesic on which the Ricci tensor 
vanishes. We therefore set a' 1 = = R^ everywhere on 
7, and for simplicity we also set aj a b to zero. 

It is known that when the Ricci tensor vanishes, the 
Riemann tensor can be decomposed in terms of a timelike 
vector u a and two symmetric-tracefree, spatial tensors 
£ a /3 and B a p (see, for example, Ref. |8|). In terms of 
frame components we have 

RaObo(u) — £ a b 

Raobc(u) = e bcd B d a , (2.64) 

Racbd(u) — 5 a b£ c d + Scd£ab — Sad£bc — 5bc£ad, 

where £ a b and B a b depend on u, are such that £ ba = £ab, 
S ab £ a b = 0, B ba = Bab, S ab B a b = 0, and e a6c is the three- 
dimensional permutation symbol. These relations can be 
substituted into Eqs. H2.28[) - l|2.30[) to g ive 

S a b = %£ab — Q a £bcQ C — ^b£ac^ C + S a b£bc^ C ^ d 

+ e acd n c B d b + EbcdSl c B d al (2.65) 

s a = £ ab n b + £ abc n b B c d n d , (2.66) 



S — £ a bQ a £l 



ar>b 



(2.67) 



In these expressions the dependence on retarded time u is 
contained in £ a b and B a b, while the angular dependence 
is encoded in the unit vector f2 a . 

It is convenient to introduce the irreducible quantities 



£* — £ a bQ a Q b , 



(2.68) 



£. 



(5 b -n n ri b )£ bc n c , 



(2.69) 



a 

£ a b = 2£ a b 2il a £bc^ 2f2;,£ ac f2 

+ {5 ab + tt a n b )£* , (2.70) 
B* a = e a b c n b B c d n d , (2.71) 



Kb = 2(s { e a -n {a n e )e b)C d^ c B c 



These are all orthogonal to fl a : £*tt a = B * a tt a = and 
£* b n b = B* ab tt b = 0. In terms of these Eqs. l|235|) - l|2~o7|) 
become 

Sab = £ab + ^a£b + £a^b + ^a^b£* 

+ B* ab + n a B* b + B* a Q b , (2.73) 
S a = £*+Q a £*+B*, (2.74) 



S = £* 



(2.75) 



When Eqs. (|2.73() - (|2.75|l are substituted into the met- 
ric tensor of Eqs. I|2.39|) ~ (|2.41|) — in which a a and u a b 
are both set equal to zero — we obtain the compact ex- 
pressions 



guu 



-1 -r 2 £* + 0(r 3 ), 



(2.76) 



Qua = -na + -r 2 {£* a +B* a )+0(r 3 ), (2.77) 



g ab = Sab-n a n b --r 2 (£; lb + B* a b) 



+ 0(r 3 ). 
The metric becomes 



2 c* 



guu = -1 - r £ 

gur 1? 



0{r 3 



{£ A + B A )+0{r% 



9AB 



r 2 n AB 



-r A {£ AB +B* AB ) 



(2.78) 

(2.79) 
(2.80) 

(2.81) 
+ 0(r 5 ) (2.82) 



after transforming to angular coordinates using the rules 
of Eq. H2.53fl . Here we have introduced the projections 



£ A = £ a Q A — £ a b^A^' , 

a, r>b _ oc n" rib 



(2.83) 

£aB = £ab^A^B = 2£ a fcf^f2 B + £*&AB, (2.84) 

B* A = B:n A = eabcn a A n b B c d n d , (2.85) 
B AB = B ab tt\n B = 2e acd n c B d b Sl a iA n b B) . (2.86) 

It may be noted that the inverse relations are £ * — £ A ^ A , 
B* a = B* A Slt £* ab = Zab^, and B* ab = B* AB ^, 
where was introduced in Eq. (|2.51() . 

The angular dependence of the quantities listed in 
Eqs. (|2.83|) - (|2.86|) can be made more explicit by express- 
ing them in terms of scalar, vectorial, and tensorial spher- 
ical harmonics. Let 



y m = {y°,Y 1c ,Y 1s ,Y 2c ,Y 2s } 



(2.87) 



be a set of real, unnormalized, spherical-harmonic func- 
tions of degree 1 = 2; explicit expressions are provided in 
Appendix A. The numerical part of the label m refers to 
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the azimuthal index m and the letter indicates whether 
the function is proportional to cos(m</>) or sm(m(f)). Vec- 
torial harmonics are defined by 



Y? = Y 



A- 



— ~ E A 



(2.88) 



where a colon indicates covariant differentiation with re- 
spect to a connection compatible with Q ab, and eab 
is the two-dimensional Levi-Civita tensor. The vectorial 
harmonics Y A have even parity, while X A have odd par- 
ity. Tensorial harmonics are defined by 

Y m fl AB , Y™ b =Y.™ b , X™ B = -X™ A . B y, (2.89) 

the harmonics Y m Q, A B and Y AB have even parity, while 
X AB have odd parity. Apart from notation and nor- 
malization, these definitions agree with those of Regge 
and Wheeler |20| , and explicit expressions appear in Ap- 
pendix A. 

We define the harmonic components S m of the tensor 
£ a b with the relations 



(2.90) 



-22 



So 


— £33 — 


Sic 


= 2S 13 , 


Sis 


= 2f 23 , 






S2c 


= \^ 


S2s 


= 2S 12 . 



Similarly, we define the harmonic components B m of the 
tensor B a b by 



B = B 33 = -(B n + B 22 ), 

B lc = 2B 13 , 

B ls = 2B 23 , 

B 2c = ^(611-622), 



(2.91) 



B-2 



2Bi 



ft is then straightforward to prove that Eqs. (|2.83(1 - (|2.86|) 

are equivalent to 

s* = ^£ m y m , 

m 

m 

£ AB = J2 £ ™( Y AB+3Y m n AB ), (2.92) 

m 

m 

B* AB = -J2^X 



AB- 



This shows that the angular dependence of these quanti- 
ties is purely quadrupolar (/ = 2). 



I. Comoving observer in a spatially-flat cosmology 

To illustrate how the formalism works we first consider 
the world line of a comoving observer in a cosmological 
spacetime with metric 

ds 2 = -dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ), (2.93) 

where a(t) is an arbitrary scale factor; for simplicity we 
take the cosmology to be spatially flat. We take the ob- 
server to be at the spatial origin of the coordinate system 
(x = y = z = 0) , and her velocity vector is given by 



u» = (1,0,0,0). 



(2.94) 



This satisfies the geodesic equation, so a^ — 0. We wish 
to transform the metric of Eq. I|2.93[) to retarded coordi- 
nates based at the world line of this observer. 

To do so we must first construct a triad of orthonormal 
spatial vectors e£. A simple choice is 

= (Co" 1 , 0,0), 



e£ = (0,0,a-\0), 
e£ = (0,0,0, a" 1 ); 



(2.95) 



these vectors are all parallel transported on 7, and we 
have LUai, — according to Eq. (|2.ip . 

Using il a — (sin 9 cos (f>, sin 9 sin (j>, cos 8) we find that 
the components of S a b, as defined by Eq. (|2.28|) . are given 
by 

Sn = —a/a + (a/ a) 2 (sin 2 9 sin 2 4> + cos 2 8), 
S12 = — (a/a) 2 sin 2 8smcj)cos(j), 
S%3 — — (a/a) 2 sin 9 cos 9 cos </), 

522 — — a/a + (a/a) 2 (sin 2 9 cos 2 cj> + cos 2 9), 

523 = — (a/a) 2 sin 9 cos 9 sin <f>, 
S 33 = — a/a + (a/a) 2 sin 2 9, 

where an overdot indicates differentiation with respect to 
t; the scale factor and its derivatives are now all evaluated 
at t = u. According to Eq. I|2.29[) . contracting S a b with 
fi b gives S a , and we obtain 

S a = -(a/a)Q a . 

Another contraction with f2 a gives 

S = —a/a, 

according to Eq. I|2.30|l . From these results it follows that 

s ab n A n b B = [-d/ a + (d/a) 2 ]n AB , 

where fl A = 3Q a /d9 A and fl A B — diag(l, sin 2 9) were 
first introduced in Sec. II G. We also have S a Vl a A = 0. 

Substituting these relations into Eqs. (|2.56|) - (|2.58|) 
shows that in the retarded coordinates, the metric com- 
ponents are given by 

= -l + r 2 (a/a) + 0(r 3 ), 



fju 



9uA = 0(r 4 ), 



(2.96) 



g AB = r 2 a^ 1 + -r 2 [a/a - (d/a) 2 ] + 0(r 3 U, 
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in addition to g ur = — 1. Not surprisingly, the metric is 
spherically symmetric. Recall that the scale factor and 
its derivatives are all functions of retarded time u. When 
the scale factor behaves as a power law, a(t) oc t a with 
a a constant, we have d/a — — a(l — a)/u 2 and d/a — 
(d/a) 2 = —a/u 2 . When instead the scale factor behaves 
as an exponential, a(t) oc e Ht with H a constant, we have 
a/a = H 2 and d/a — (a/a) 2 = 0. 

To help clarify the meaning of these results, we present 
next an ab initio derivation of Eq. I|2.96[) . We take the 
metric of Eq. (|2.93|) and switch to conformal time 77, 
which is defined by the relation drj — dt/a(t). The metric 
becomes 



ds 2 



a 2 (rj)(-drf + dx 2 + dy 2 + dz 2 ). 



We then introduce spherical coordinates (p, 9, <p) through 
the relations x = p sin 9 cos </>, y = p sin 9 sin if), z = 
pcos9, and the null coordinate u — r\ — p. The metric 
now reads 



ds 2 



a 2 (u + p)(-du 2 - 2 dudp + p 2 dn 2 ) , (2.97) 



where dfl 2 = Q AB d9 A d9 B = d9 2 + sin 2 9 d<j) 2 . While u is 
a retarded-time coordinate and p is a radial coordinate, 
these are distinct from u and r, and Eq. (|2.97() does not 
match the form of Eq. (|2.96|l . Since u and u can both be 
used to label light cones centered at p = 0, there must 
exist between them a relation of the form u = u(u). And 
since r and p can both be used to parameterize the null 
generators of a light cone u — constant (although p is 
not an affine parameter), there must exist between them 
a relation of the form r = r(u,p). We shall now obtain 
these relations. We recall that the world line is located 
at p = r = 0. 

In the coordinates (u,p,9,(f>), the observer's velocity 
vector is given by 



1 



a(u) 



,0,0,0 



(2.98) 



where the scale factor is evaluated at p = and expressed 
as a function of u only. The null vector k a = —d a u is tan- 
gent to the null cones u — constant, and its components 
are given by k a = (0, a~ 2 , 0, 0), where the scale factor is 
now a function of u + p. We have k^u^ = —l/a(u). 

We are seeking a null coordinate u and a null vector 
field k a = —d a u such that k^u^ — —1; refer back to 
Eqs. H2.19JI and l|2.25|l . It is easy to see that this should 
be given by k a = a(u)k a . We therefore define u with the 
statement 



and 



du = a(u) du, 



a 2 (u + p) ' 



(2.99) 



(2.100) 



is tangent to the light cones u = constant. This vector 
satisfies the geodesic equation k"^^ = 0, and the affine 



parameter on the null generators is r. From Eq. (|2.100|) 
we have dp/dr — a(u)/a 2 (u + p), which integrates to 



p a 2 (u + p') 
a(u) 



dp', 



(2.101) 



taking into account the boundary values r(u, p = 0) = 0. 
Equations l(!H)9) and l|2.101|) give us the transformation 
between the old coordinates (u, p) and the new coordi- 
nates (u, r). 

After applying the coordinate transformation to the 
metric we obtain 



ds 2 



v a 2 {u + p) 2r, D 



a 2 (u) a(u) 



du 2 -2 dudr+a 2 (u+p)p 2 dtt 2 , 



(2.102) 

where r u = dr/du. To show that this matches the results 
of Eq. (|2.96|l , we must evaluate the integral of Eq. I|2.101|l . 
It is sufficient to work in a neighborhood of p = 0, and 
a 2 (u + p) can be expressed as a Taylor expansion. This 
yields 



ap 



1 



a 2 " 



0(P 3 



where a prime indicates differentiation with respect to u; 
here and below, the scale factor and its derivatives are 
evaluated at p — and expressed as functions of u only. 
From this we gather that r a = a'p + a"p 2 + 0(p 3 ), and 
Eq. H2.102I) gives 



and 
9ab 



a 2 p 2 n 



2-p- 
a 



p 2 + 0(p*) 



a 2 " 



0{pf 



Expressing these results in terms of r instead of p, 
and converting ^-derivatives into it-derivatives using 
Eq. (|2~g9"|) . returns the resu lts of Eq. lf23rJ|l . It should be 
noted that while Eq. $H)^i gives the metric in a neigh- 
borhood of r = 0, the expression given in Eq. (|2.102|) 
holds globally. 



J. Circular motion in Schwarzschild spacetime 

As another example we consider the world line of 
a freely-moving observer in circular motion around a 
Schwarzschild black hole. In the usual Schwarzschild co- 
ordinates (t s , r s ,9 Sl 4> s ) the metric is given by 

ds 2 = -(1 - 2M/r s ) dt 2 s + (1 - 2M/r s )- 1 dr 2 s 

(2.103) 



r 2 (d9 2 + sin 2 I 



*2). 



where M is the mass of the black hole. The observer 
moves on a circular orbit of radius r s = R with an angular 
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velocity d(p s /dt s = f2 = yj M / R 3 , in the equatorial plane 
9 S = 7r/2. The velocity vector is 



u" = 7 (i,o,o,n), 



(2.104) 



where 7 = (1 — 3A//i?) -1 / 2 is a normalization factor. The 
motion is geodesic, and we can once more set a M = 0. 
Because Rap = for the Schwarzschild spacetime, we 
will rely on the results presented in Sec. II H. 
The vectors 



(0, ft 0,0), 



£ = (0,0,l/i?,0), 



and 



s£ = 7(0^7/3, 0,0,ftAR), 



where (3 — (1 — 2M/i?) 1 / 2 , are normalized, mutually or- 
thogonal, and all orthogonal to u M . As such they form a 
valid set of spatial vectors, but this choice is not op- 
timal because except for eg, the vectors are not par- 
allel transported on the world line. By forming linear 
superpositions, however, and choosing the coefficients 
appropriately, we can find a set of parallel-transported 
vectors e^. We choose = cos(f2-r)e{? — sin(f2i-)e|f, 



sin(rir)e^ + cos(f2r)e^, and e 



M _ 
3 — 



or 



jttR fty \ 

— sin$,/5cos$,0,-— sin$ , 
P R J 

{^^- cos $, sin 0, ^ cos , (2.H M i 



R 



1 



0,0,--, 



where f3 = (1 - 2M/R) 1 ' 2 , j = (1 - 3M/R)- 1 / 2 , and 

$ = ftr. (2.106) 

The vectors are all parallel transported on the world 
line, and according to Eq. (|2.1|) . we have cj a ;, = 0. 

The electric part £ a b of the Riemann tensor is defined 
by Eq. (|2.64ll , and Eq. I|2.90(l gives its decomposition into 
harmonic components £ m . Using the tetrad introduced 
previously we find that the nonvanishing components are 



So 



M 



£lr. — — 



R 2 (R-3M)' 
3M R — 2M 



£■2 



2R 3 R 
3M R 



3M 
2M 



■ cos 2$, 



(2.107) 



2R 3 R - 3M 



■ sin 2$. 



The magnetic part B a b of the Riemann tensor is defined 
in Eq. 12.64J1 and decomposed in harmonic components 
in Eq. (j2.91|) . Its nonvanishing components are 



6MCI 



, 2M ^ 

. , 1 COS 9, 

R{R - 3M) V R 



(2.108) 



6MCI 



2M 

R(R^my 1 '^ s ' m " > - 



Equivalent results were obtained by Alvi |2lJ . based on 
earlier work by Fishbone and Marck |23j ] . 

To construct the metric we must form the quantities 



£*, £* A , £* AB , B* A , and B* AB , defined by Eqs. (|2.68[l and 
(|2~k^ - (|2~fSjl . For this we use Eq. $2~§2^ and the spherical 
harmonics listed in Appendix A. We obtain 



£* 



M 



2R 2 {R-3M) 



3cos^(9 - 1) 



£t = - 



3MR-2M . 2n . 

— - sin ft cos 2 

2R 3 R — 3M v 

3M \R 



2R 3 (R-3M) 



(R - 2M) cos 2(0 - $)] sin 9 cos 9, 



£1 = 



3M R-2M 
2R 3 R - 3M 

3M 
2R 2 (R-3M) 



sin 2 6* sin 2(0- $), 



sin 2 9 



3M R - 2M 



£94 



2R 3 R - 3M 
3M R - 2M 



-(l + cos 2 9) cos 2(0- $), 



R 3 R-3M 



sin 9 cos 9 sin 2(0 — $), 



3M 



2R 2 {R-3M) 
3M R-2M 
h 2R 3 R - 3M 



sin 2 6»d + cos 2 s ) cos2( 



(2.109) 



and 
K 

Bmj, 



3Mtt 



2M „ , 
i?(i?-3A/)V 1 "lT COS0Sm( 



3MO / 2M 
i?(i? - 3 M) V ~ 
x sin 6»(l - 2 cos 2 0) cos(0 - $), 



6MO 



,,1 sin W sin 

R{R - 3M) V i? V 



6MO / 2M . 2 



6Mf2 



i?(i? - 3M) 



1 



2M 



sin t/ sin 



(2.110) 



In these expressions, the components of the Riemann ten- 
sor are all evaluated at r = u, so that <& = £lu. Substitut- 
ing them into Eqs. H2.79[H2.82[) gives the Schwarzschild 
metric in the retarded coordinates (it, r, 0,0); these are 
based at the world line of an observer moving on a cir- 
cular orbit of radius R with an angular velocity Q = 

yfWJW. 
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III. MOTION OF A SMALL BLACK HOLE IN 
AN EXTERNAL UNIVERSE 

A. Matched asymptotic expansions 

In this section we consider a nonrotating black hole of 
small mass m moving in a background spacetime with 
metric g a /3, and we seek to determine the equations that 
govern its motion. We will employ the powerful technique 
of matched asymptotic expansions 0, 0, 0, 0, Q, 0, 
EH ] and make use of the retarded coordinates developed 
in Sec. II. 

The problem presents itself with a clean separation of 
length scales, and the method relies heavily on this. On 
the one hand we have the length scale associated with 
the small black hole, which is set by its mass m. On 
the other hand we have the length scale associated with 
the background spacetime, which is set by the radius of 
curvature TZ; this is defined so that a typical component 
of the background spacetime's Riemann tensor is equal 
to 1/1Z 2 up to a numerical factor of order unity. We 
demand that m/lZ <C 1. For simplicity we assume that 
the background spacetime contains no matter, so that 
its metric is a solution to the Einstein field equations in 
vacuum. 

Let r be a meaningful measure of distance from the 
small black hole, and let us consider a region of spacetime 
defined by r < 7-j, where is a constant that is much 
smaller than TZ. This inequality defines a narrow world 
tube that surrounds the black hole, and we shall call 
this region the internal zone. In the internal zone the 
gravitational field is dominated by the black hole, and 
the metric can be expressed as 

g(internal zone) = g (black hole) + H 1 /K + H 2 /1Z 2 H , 

where g (black hole) is the metric of a nonrotating 
black hole in isolation (as given by the unperturbed 
Schwarzschild solution), while Hi and H 2 are corrections 
associated with the conditions in the external universe. 
The metric of Eq. I|3.1|l represents a black hole that is 
distorted by the tidal gravitational field of the external 
universe, and Hi, H 2 are functions of the spacetime co- 
ordinates that can be obtained by solving the Einstein 
field equations. They must be such that the spacetime 
possesses a regular event horizon near r = 2m, and such 
that g (internal zone) agrees with the metric of the ex- 
ternal universe — the metric of the background space- 
time in the absence of a black hole — when r » ra. As 
we shall see in Sec. Ill B, Hi actually vanishes and the 
small correction H 2 /TZ 2 can be obtained by employing 
the well-develo ped tools of black-hole perturbation the- 
ory EE0MII1. 

Consider now a region of spacetime defined by r > r e , 
where r e is a constant that is much larger than m; this re- 
gion will be called the external zone. In the external zone 
the gravitational field is dominated by the conditions in 



the external universe, and the metric can be expressed as 

g (external zone) = g (background spacetime) 

+ mh 1 +m 2 h 2 -) , (3.2) 

where g (background spacetime) is the unperturbed met- 
ric of the background spacetime in which the black hole is 
moving, while hi and hi are corrections associated with 
the hole's presence; these are functions of the spacetime 
coordinates that can be obtained by solving the Einstein 
field equations. We shall truncate Eq. (|3.2|) to first or- 
der in m, and mhi will be calculated in Sec. Ill C by 
linearizing the field equations about the metric of the 
background spacetime. 

The metric g (external zone) is a functional of a world 
line 7 that represents the motion of the small black hole 
in the background spacetime. Our goal is to obtain a 
description of this world line, in the form of equations of 
motion to be satisfied by the black hole; these equations 
will be formulated in the background spacetime. It is im- 
portant to understand that fundamentally, 7 exists only 
as an external-zone construct: It is only in the external 
zone that the black hole can be thought of as moving on 
a world line; in the internal zone the black hole is re- 
vealed as an extended object and the notion of a world 
line describing its motion is no longer meaningful. 

Equations jSJJ and (|3.2|l give two different expressions 
for the metric of the same spacetime; the first is valid in 
the internal zone r < 7"j <C TZ, while the second is valid 
in the external zone r > r e 3> m. The fact that 1Z m 
allows us to define a buffer zone in which r is restricted to 
the interval r e < r < T{ . In the buffer zone r is simultane- 
ously much larger than m and much smaller than 1Z — a 
typical value might be VmTZ — and Eqs. |3. Ij3.2|l are 
simultaneously valid. Since the two metrics are the same 
up to a diffeomorphism, these expressions must agree. 
And since g (external zone) is a functional of a world line 
7 while g(internal zone) contains no such information, 
matching the metrics necessarily determines the motion 
of the small black hole in the background spacetime. 

Matching the metrics of Eqs. (|3.1|) and (|3.2|) in the 
buffer zone can be carried out in practice only after per- 
forming a transformation from the external coordinates 
used to express g (external zone) to the internal coordi- 
nates employed for g(internal zone). The details of this 
coordinate transformation are presented in Appendix C 
and the end result of matching is revealed in Sec. Ill D. 

B. Metric in the internal zone 

To proceed with the program outlined in the previous 
subsection we first calculate the internal-zone metric and 
replace Eq. (|3.1|) by a more concrete expression. We re- 
call that the internal zone is defined by r < r.; < R, 
where r is a suitable measure of distance from the black 
hole. 

We begin by expressing g (black hole), the 
Schwarzschild metric of an isolated black hole of 
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mass m, in terms of retarded Eddington-Finkelstein 
coordinates (u,f, 9 A ), where u is retarded time, f the 
usual areal radius, and 9 A — (9, </>) are two angles on the 
spheres of constant u and f. The metric is given by 



ds 2 



-fdu 2 -2dudf + f 2 dn 2 



where 



/ = 1 



2m 
r 



sin 



(3.3) 



(3.4) 



is the line 



and dVL 2 = n AB d9 A d9 B = d9 2 
element on the unit two-sphere. In the limit r»m this 
metric achieves the asymptotic values g uu — ► — 1, g a f = 
— 1) 9uA = 0; an d 9 AS = ^ 2 &AB', these are appropriate 
for a black hole immersed in a flat spacetime charted by 
retarded coordinates. 

The corrections Hi and H 2 in Eq. 1|3.1|) encode the 
information that our black hole is not isolated but in 
fact immersed in an external universe whose metric be- 
comes g (background spacetime) asymptotically. In the 
internal zone the metric of the background spacetime can 
be expanded in powers of f/TZ and expressed in a form 
that can be directly imported from Sec. II. If we assume 
that the "world line" f = has no acceleration in the 
background spacetime (a statement that will be justified 
shortly), then the asymptotic values of g (internal zone) 
must be given by Eqs. (ZZSHZBZ) : 

gnu -l-r 2 £* + 0(f 3 /TZ 3 ), 

guf 1) 



guA 

gAB 
where 



~f 3 (£* A + B* A ) +0{f A /n% 

f 2 (l AB - \f\£\ B + B* AB ) + 0(f 5 /K 3 ), 



£* 

t-A 

c* 
C AB 



£ab& a & , 
£ab& A & j 

2£ a b& A & B + £*&AB 



and 



B A = £ a b c O, A Q B c d Q , 
B* AB = 2e acd n c B d M0 



are the tidal gravitational fields that were first introduced 
in Sec. II H. Recall that fl a = (sin 9 cos (f>, sin 9 sin </>, cos 9) 
and Q, A = dVL a /d9 A . Apart from an angular depen- 
dence made explicit by these relations, the tidal fields de- 
pend on u through the frame components £ a b = RaObo = 
0(1/TZ 2 ) and B\ = \e acd R obcd = 0(1/TZ 2 ) of the Rie- 
mann tensor. (This is the Riemann tensor of the back- 
ground spacetime evaluated at f = 0.) Notice that we 
have incorporated the fact that the Ricci tensor vanishes 
at f = 0: the black hole moves in a vacuum spacetime. 



The modified asymptotic values lead us to the follow- 
ing ansatz for the internal-zone metric: 

gun = -f[l + f 2 ei (f)£*]+0(f 3 /TZ 3 ), (3.5) 
gar = -1, (3.6) 

guA = f 3 [e 2 (f)£ A + b 2 (f)B A ]+0(f 4 /K%(3.7) 



g AB = f 2 Cl AB - -f 4 [e 3 (f)£* AB +b 3 (f)B AB ] 
+ 0(f 5 /TZ 3 ). 



(3.8) 



The five functions e\, e 2 , b 2 , and b 3 can all be deter- 
mined by solving the Einstein field equations; they must 
approach unity when r 3> m and be well-behaved at 
r = 2m (so that the tidally distorted black hole will have 
a nonsingular event horizon). In Appendix B, I show that 
they are given by 



ei(r) = e 2 (r) = /, 



e 3 (r) 



1 



2m 2 



and 



62 (r) = /, 63(f) = 1. 



(3.9) 



(3.10) 



It is clear from Eqs. H3.5f) - H3.8f) that the assumed de- 
viation of g (internal zone) with respect to g (black hole) 
scales as 1/1Z 2 . It is therefore of the form of Eq. H3.1(l 
with Hi = 0. The fact that Hi vanishes comes as a conse- 
quence of our previous assumption that the "world line" 
f = has a zero acceleration in the background space- 
time; a nonzero acceleration of order 1/1Z would bring 
terms of order 1/7?. to the metric, and Hi would then 
be nonzero. The perturbed metric of Eqs. (|3.5I) - (|3.10() 
differs from the one presented by Detweiler |15| only 
by a transformation from Schwarzschild to Eddington- 
Finkelstein coordinates, and a transformation from the 
Zerilli gauge |24| gauge adopted by him to the retarded 
gauge adopted here. 

Why is the assumption of no acceleration justified? As 
I shall explain more fully in the next paragraph, the rea- 
son is simply that it reflects a choice of coordinate sys- 
tem: setting the acceleration to zero amounts to adopting 
a specific — and convenient — gauge condition. 

Inspection of Eqs. I|3.5|l l|3.8|l reveals that the angu- 
lar dependence of the metric perturbation is generated 
entirely by scalar, vectorial, and tensorial spherical har- 
monics of degree 1 — 2; this observation was elaborated 
toward the end of Sec. II H and in Appendix A. In par- 
ticular, H 2 contains no I = and / = 1 modes, and this 
statement reflects a choice of gauge condition. Zerilli 
has shown |24| that a perturbation of the Schwarzschild 
spacetime with I = corresponds to a shift in the mass 
parameter. As Thorne and Hartle have shown |8j, a 
black hole interacting with its environment will undergo 
a change of mass, but this effect is of order m 3 /TZ 2 and 
thus beyond the level of accuracy of our calculations. 
There is therefore no need to include I = terms in H 2 . 
Similarly, it was shown by Zerilli that odd-parity per- 
turbations of degree I = 1 correspond to a shift in the 
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black hole's angular-momentum parameters. As Tliorne 
and Hartle have shown, a change of angular momentum 
is quadratic in the hole's angular momentum, and we can 
ignore this effect when dealing with a nonrotating black 
hole. There is therefore no need to include odd-parity, 
I = 1 terms in H2. Finally, Zerilli has shown that in a 
vacuum spacetime, even-parity perturbations of degree 
I = 1 correspond to a change of coordinate system — 
these modes are pure gauge. Since we have the freedom 
to adopt any gauge condition, we can exclude even-parity, 
I = 1 terms from the perturbed metric. This leads us 
to Eqs. H3.5fl - H3.8fl . which contain only I = 2 perturba- 
tion modes; the even-parity modes are contained in those 
terms that involve £ a b, while the odd-parity modes are 
associated with B a b- The perturbed metric contains also 
higher multipoles, but those come at a higher order in 
1/7?.; for example, the terms of order 1/7\L 3 include / = 3 
modes. We conclude that Eqs. (|3.5II - (|3.8I) is a sufficiently 
general ansatz for the metric in the internal zone. 

It shall prove convenient to transform g (internal zone) 
from the quasi-spherical coordinates (f, 9 A ) to a set of 
quasi-Cartesian coordinates x a = rti a {Q A ); the transfor- 
mation rules are worked out in Sec. II G. This gives 

gun = -f{l + f 2 f£*)+0(f 3 /K 3 ), (3.11) 

gu & = -n a + ^f 2 f(£* a + B* a ) + 0(f 3 /K 3 ), (3.12) 

1 / 777 2 \ - 1 

g- ab = 5 ab -n a n b --f 2 \l-2-^\£* ab --r 2 B* ab 

+ 0(f :3 /n 3 ), (3.13) 
where / = 1 — 2m /f and where the tidal fields 
£* = £ a b& a & b , 

£* a = (s a b ~n a n")£ bc n c , 

&ab = 2£ a b — 2Vl a £ bc VL c — 2fl b £ ac il C + (S ab + V, a il b )£* 7 

B* a = e abc n b B c d n d , 

B* ab = s acd n c B d e (s e b - n e n b ) + e bcd n c B d e (s e a - n e n a ), 

were first introduced in Sec. II H. The metric of 
Eqs. H3.11|l - H3.13|) represents the spacetime geometry of 
a black hole immersed in an external universe and dis- 
torted by its tidal gravitational field. 

C. Metric in the external zone 

We next move on to the external zone and seek to re- 
place Eq. i|3.2[l by a more concrete expression; recall that 
the external zone is defined by m < r e < r. We take 
advantage of the fact that in the external zone, the grav- 
itational perturbation associated with the presence of a 
black hole cannot be distinguished from the perturbation 
produced by a point particle of the same mass. 

The external-zone metric is decomposed as 

ga/3 = fJa/3 + h af3 , (3-14) 



where g a fj is the metric of the background spacetime and 
h a fj — 0(m) is the perturbation; we shall work consis- 
tently to first order in m and systematically discard all 
terms of higher order. We relate h a p to trace-reversed 
potentials y a p, 

Kp = 7 Q/3 - - (fl^Ty*)^' ( 3 - 15 ) 

and we impose the Lorenz gauge condition 

7°*% = 0; (3.16) 

indices are raised and lowered with g a @ and g a p, re- 
spectively. With the understanding that the background 
spacetime contains no matter, linearizing the Einstein 
field equations produces the wave equation 

□ 7 Q/3 + 2i? Q y /3 5 7 7 ' 5 = -16irT af3 (3.17) 

for the potentials. Here, □ = g a/3 V ' Q V ' p is the wave 
operator and 

T a/3 (x)=m j g a ll {x,z)g l3 v {x,z)u tl u"5 i {x,z)dT (3.18) 

is the stress-energy tensor of a point particle of mass 
m traveling on a world line 7; 8i{x,z) is a scalarized, 
four-dimensional Dirac functional, the world line is de- 
scribed by relations z^(t) in which t is proper time, and 
= dz^/dr is the particle's velocity vector. Solving the 
linearized field equations produces 

7a/3(a;)=4m / G a p ltv (x,z)u ,i u u dr, (3.19) 
J 1 

where G a /3 tll/ (x, z) is the retarded Green's function 2 
associated with Eq. I|3.17|l . 

We now place ourselves in the buffer zone (where m <C 
r«R and where the matching will take place) and work 
toward expressing g (external zone) as an expansion in 
powers of r/TZ. For this purpose we adopt the retarded 
coordinates (u, r£l a ) of Sec. II and rely on the machinery 
developed there. 

We begin with g a p, the metric of the background 
spacetime. We have seen in Sec. II H that if the world 
line 7 is a geodesic, if the vectors are parallel trans- 
ported on the world line, and if the Ricci tensor vanishes 
on 7, then the metric takes the form given by Eqs. l|2~7o]l- 
(|2.78l) . This form, however, is too restrictive for our pur- 
poses: We must allow 7 to have an acceleration, and allow 
the basis vectors to be transported in the most general 



2 The normalization of the gravitational Green's function varies 
from author to author. Here the normalization is such that the 
Green's function obeys a wave equation with a right-hand side 
given by — 47r</^ (x, z)g^} (x, z)S.i(x, z). The factor of 4ir accounts 
for the factor of 4 on the right-hand side of Eq. 13.191 . 
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where r< and r> are the values of the proper-time pa- 
rameter at which 7 enters and leaves jV(x) , respectively. 
The third integration contributes nothing because x is 
then in the past of z{f) and the retarded Green's func- 
tion vanishes. For the second integration, x is the normal 
convex neighborhood of z(t), and the retarded Green's 
function can be put in the Hadamard form 0, 0, I27I |2S| 



u 



FIG. 2: The region within the dashed boundary represents 
the normal convex neighborhood of the point x. The world 
line 7 enters the neighborhood at proper time t< and exits at 
proper time r>. Also shown is the retarded point x' = z(u) 
and the null geodesies that links x to the world line. 



G a pnu(x, z) = U a p^{x, z)S + {a) + V aj3a y{x, z)6 + (-a), 

(3.23) 

where U a /3 flI/ (x, z) and V a (}^ v {x, z) are smooth bitensors 
3 , a(x, z) is Synge's world function 0,0], and 8 + {<r) and 
0+(c) are Dirac and Heaviside distributions restricted 
to the future of z(t) 0, l28| . To integrate over the 
Dirac term we change variables from r to a, noticing 
that a increases as z(t) passes through the retarded 
point z(u) = x'; recalling the definition of retarded 
distance given in Eq. I|2.1U[1 . the integral evaluates to 
Uapj'S'iXjX 1 )^ u s jr. The integration over the Heavi- 
side term is cut off at r = u, and we obtain our final 
expression for the perturbation: 



4m 



7«/?(z) = — U aPl , s ,{x,x')if V +j t a f(x). (3.24) 



way compatible with their orthonormality property; this 
transport law is given by Eq. I|2.1jl . 



dr 



a a u^ + uj a e£, 



where a a (r) = a^e^ are the frame components of the 
acceleration vector = Du^ /dr, and w a b(r) = — uibaij) 
is a rotation tensor to be determined. Anticipating that 
a a and oj a b will both be proportional to to, we express 
the metric of the background spacctime as 



Qua 



9ab 



-1 - 



2ra a n a ~r 2 £* + 0(r 3 /lZ 3 ), (3.20) 
f r(6 a b -n a n b )a b -ru; ab n b 

3 /<r»3-| 



+ ^r 2 (£:+B:)+0(ryn% 



Sab 



ila^b 
„3 /<T>3\ 



1 



'iCb + K 



+ 0(r 3 /TZ 3 ), 



(3.21) 



(3.22) 



where £*, £*, £* b , 0*, and B* ab are the tidal gravita- 
tional fields first introduced in Sec. II H. The metric of 
Eqs. (|3.20|) - (|3.22l) is obtained from the general form of 
Eqs. (|2.39|l - (|2.41|) by neglecting terms quadratic in a a 
and uj a b and specializing to a zero Ricci tensor. 

To express the perturbation h a p as an expansion in 
powers of r/TZ we assume that x is sufficiently close to 
7 that a portion of the world line traverses Af(x), the 
normal convex neighborhood of the point x (see Fig. 2); 
this assumption is compatible with the condition r<K. 
We then re-express Eq. (|3.19() as 



7a/3 



4m 



GapuvU^u 1 ' dr, 



Here, primed indices refer to the retarded point x' = z(u) 
associated with x, and 



lit( x ) = 4m / V a p ia ,(x,z)vru v drr 

+ 4:m / G a fj^(x, z)u^v, v dr 



= 4m / G a0fiu (x, z)u»u v dr (3.25) 
J —00 

is the "tail part" of the gravitational potentials. No- 
tice that I have introduced a short-hand notation in the 
last line of Eq. I|3.25ll : the important point is that the 
integral of G a ff fjtu (x, z^u^v^ up to r = u~ avoids the sin- 
gular behavior of the Green's function on the light cone 
a(x, z) = 0. 

We must next express Eq. I|3.24|l in the form of an 
expansion in powers of r/TZ. For this we shall need the 
expansion 4 



U^ys'U^'u 5 ' = g%9 P ' p) \u a >u p , + 0{r 3 /TZ 3 ) 



which contains no terms at order r/TZ or (r/TZ) 2 . We 
shall also need the general expansion of a tensor A a p (x) 



3 The tail part of the Green's function is denoted V a ^ M ^(x, z) 
in this work, but most authors insert a minus sign and call it 
-v a pft V (.x,z). 

4 This result is implicitly contained in App end ix A of Ref . . It 
is derived from scratch in Sec. 15 of Ref. Il7l . 
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/ 5. 



in terms of its values at a neighboring point x 



A a p{x) = g°ig [ A a ,p, - Aa>p>. n , a r + 



Here A a p will stand for y^ir and x' will be the retarded 
point z(u) associated with x; accordingly, Eq. i|2.11[l gives 
a a = —r(u a + tt a e" ). Combining all these results, 
Eq. H3.24jl becomes 



4m 



+ 0(mr 2 /K 3 ) 



(3.26) 



where 7^g/ is the tensor of Eq. 1|3.25[) evaluated at x', 
and 



7 ^ Y (x') = 4m / V y G^- F (i',z)//rfT (3.27) 
emerges during the computation of 



Vy7«73' = 4m{d Y u)V a ,p, 5 , t ,u d u e +j^, Y ; 

the term proportional to dy>u disappears after contrac- 
tion with a 1 . 

At this stage we introduce the fields 



h% = 4m 



x vFu" dr, 



(3.28) 



r tail 
' l a'/3'7' 



4m 



x u^u v dr 



(3.29) 



and recognize that the metric perturbation obtained from 
Eqs. (|3~T5) and fiTZty is 



Kp{x) = g a a g p f3 



2m f 

\2u a 

r V 

+ (9(mr 2 /^ 3 ) 



U/3' + .9a'/3' ) + h^p, 



(3.30) 



This is the desired expansion of the metric perturbation 
in powers of r/TZ. Our next task will be to calculate the 
components of this tensor in the retarded coordinates 
{u,rfl a ). 



5 This is Eq. (A21) of Ref. This result is derived from scratch 
in Sec. 5 of Ref. I??! . 



The first step of this computation is to decompose h a p 
in the tetrad (eg,e") that is obtained by parallel trans- 
port of (u a , e" ) on the null geodesic that links x to its 
corresponding retarded point x' = z(u) on the world line. 
(The vectors are parallel transported in the background 
spacetime.) The projections are 



h 00 (u,r,n a ) 



h ob {u,r, n a ) 



h ab {u,r, Q a ) 



h a0 e%e% = ^ + hlf(u) 

+ r[h^(u) + hlf c (u)n c ] 
+ 0(mr 2 /K 3 ), 
Kpe^ = hlf(u) 
+ r[h^(u) + hlt(u)n c ] 
+ 0(mr 2 /K 3 ), 



(3.31) 



(3.32) 



2m 



-6 ab + h a f(u) 



+ r[h a ^(u) + h a f c (u)n c } 

+ 0(mr 2 /K 3 ). (3.33) 

On the right-hand side we have the frame components 
of h^p, and hffip, , taken with respect to the tetrad 

(u a ' , e"' ) ; these are functions of retarded time u only. 
The perturbation is now expressed as 

h a ff = /iooe°e° + ft, b(e°e^ + e b a e° ) + h a be^e b 

and its components are obtained by involving Eqs. (|2.37() 
and (|2.38|) . which list the components of the tetrad vec- 
tors in the retarded coordinates; this is the second (and 
longest) step of the computation. Noting that a a and 
uj a b can both be set equal to zero in these equations (be- 
cause they would produce negligible terms of order m 2 
in hap), and that S ab , S a , and S can all be expressed in 
terms of the tidal fields £*, £*, £* b , B* a , and B* ab using 
Eqs. (|2~73ll - lj2~75|l . we arrive at 



2m 
r 



h 



tail 
00 



(2m£* 



■tail 
'000 



0(mr 2 /TZ 3 



(3.34) 



r 



h 



tail 
0a 



2m£*n a 



+ —(£* 



■B* 



l tail , q 7 tail , i tailQfc 
"OaO + "a^-OOO + "Oah" 



O /7 tall O 
"a'^OOb" 



+ 0(mr 2 /K 3 ), 



(3.35) 



2m 



(S a b + ^a^b) + ^a^b^oO 1 ' + ^a^Qh + ^ft^Oa 



'tail 



n a £* b +£:n b + B* a 



n a B* b + a b B a ) +Q a n b (h [ 
n a (h 



tail 
000 



tail 
L 0b0 



^tail 
30c 

tail 



i tailQC 
"•OOc" 



(h u 



tail 
bO 



+ 0{mr 2 /K 3 ). 



(3.36) 
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These are the coordinate components of the metric per- 
turbation h a fj in the retarded coordinates (u,rQ a ), ex- 
pressed in terms of frame components of the tail fields 
tfa^pi an d tf^p, The perturbation is expanded in pow- 
ers of r/TZ and it also involves the tidal gravitational 
fields of the background spacetime. 

The external-zone metric is obtained by adding g a p 
as given by Eqs. l|3.2L)[l - il3.22[l to h a p as given by 
Eqs. (E33-(E311l- The final result is 

gnu = -l-r 2 £* +0(r 3 /1l 3 ) 
2m 



tail 
00 



(2m£* - 2a a tt a + h, 



tail 
(]()() 



+ Ca^)+0(mr 2 /^ 3 ), 
gua = -n a + ^r 2 (£:+B:)+0(r 3 /n 3 ) 



(3.37) 



2m. 



O 4- /? tail -I- O h 



tail 
00 



2m£*n a 



2m 

—(s* a + B:)+(s a b - n a n b )a b - cu ab n b 



, Ltail , q ttail , LtailQb , q itailQb 

+ 0(mr 2 /K 3 ), 

■ib = S a b 



(3.38) 



1 



n a n b --r'{£* ab + B* ab )+0(r d /n 6 ) 



2m 



+ — {Sab + n a n b ) + n a ttbhm + fl a h 



,tail 
Ob 



+ Q b h t f + h t af + r 



— [£* +n r 

2 \ c ab + ! ' L aO b 



+ £*n b + B* ab + n a B* b + n b B* a j + fi„Q b (/^d 



( I, tail i ^tail 
\ n ab0 



itailQC 



+ 0{mr 2 /K 3 ). 



(3.39) 



D. 



Matching: motion of the black hole in the 
background spacetime 



Com parison of Eqs. j!TTT|) - (|3~T3|) and Eqs. ^H7|> - 
(|3.39l) reveals that the internal-zone and external-zone 
metrics do no match in the buffer zone. But as the 
metrics are expressed in different coordinate systems, 
this mismatch is hardly surprising. A meaningful com- 
parison of the two metrics must therefore come after a 
transformation from the external coordinates (u,r£l a ) to 
the internal coordinates (u, fCl a ). This transformation is 
worked out in Appendix C, and it puts the external-zone 
metric in its final form 



-1 - f 2 £ * + 0(f 3 /K 3 ) 

2™ - r . ^* 

+ — + r 4m£* - 
r 

+ 0{mf 2 /Tl 3 ), 



1 ttail 
2 ,l 00a 



(3.40) 



w-a = -n a + -f 2 (£* a 



4m 



1 t tail 
^"■OOb 



(4* 



-K 

b*: 



- 0(f 3 /K 3 ) 

{d a b - n a n b ) [ a b 



tail 
'0b0 



(UJab 



j tail \qI) 



+ 0(mf 2 /Tl 3 ), 

lab = Sab - & a n b - 



(S* ab + B* ab ) 
+ 0(f 3 /K 3 ) + 0{mf 2 /K 3 ). 



(3.41) 



(3.42) 



Except for the terms involving a a and u) ab , this metric is 
equal to g(internal zone) as given by Eqs. (|3.11|K(|3.13|) 
linearized with respect to m. 

A precise match between g (external zone) and 
g (internal zone) is produced when we impose the rela- 
tions 



_ 1 l, tail 
— 2 00a 



'tail 
*0a0 



and 



While Eq. 



, , i tail 

U ab — ft [a&] 



(3.43) 



(3.44) 



tells us how the black hole moves in 
the background spacetime, Eq. I|3.44ll indicates that the 
vectors are not Fermi- Walker transported on the world 
line. 

The black hole's acceleration vector a p = a a e a l can 
be constructed from the frame components listed in 
Eq. H3.43fl . A straightforward computation gives 

^ = + (2hlf p - ft$£)uV, (3.45) 



where the tail integral 



itail 



4m 



x vP'u"' ' dr' 



(3.46) 



was previously defined by Eq. (|3.29l) . Here, the unprimcd 
indices refer to the current position z{j) on the world line, 
while the primed indices refer to a prior position z(t ); 
the integral is cut short at r' = t~ in the manner defined 
by Eq. I|3.25(l . These are the MiSaTaQuWa equations of 
motion, as they were first presented by Mino, Sasaki, and 
Tanaka 01 > an d later rederived by Quinn and Wald [ill • 
Substituting Eqs. I|3.43[) and (|3.44|) into Eq. (|2.1|) gives 
the following transport equation for the tetrad vectors: 



Dett 



dr 



l^(2hif p 



ti?R)u v e x u p 



v vp\j 

+ (g»P+u>*u»)h^ p] u»e x 



dr 



(3.47) 

\ (u^u" + ff x eP a - gTefy'tift, (3.48) 



This can also be written in the alternative form 
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that was first proposed by Mino, Sasaki, and Tanaka. 
Both equations state that in the background spacetime, 
the tetrad vectors are not Fermi- Walker transported on 
7; the rotation tensor is nonzero and given by Eq. 1)3. 44f) . 
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APPENDIX A: SPHERICAL HARMONICS 

In this Appendix I provide explicit expressions for the 
scalar, vectorial, and tensorial spherical harmonics intro- 
duced in Sec. II H. All harmonics are of degree I = 2. 

The scalar harmonics are 



Y° 




COS 2 1 
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Y lc 


— sin 6 


'cos t 


' COS 0, 


Y ls 


= sin 6 


] cost 


?sin0, 


Y 2c 


= sin 2 


9 cos 


20, 


Y 2s 


= sin 2 


9 sin 


2(f>. 



The vectorial harmonics are defined by YT = X™ (even 
parity) and X™ = —e^Y^ (odd parity), where a colon 
indicates covariant differentiation with respect to a con- 
nection compatible with Qab = diag(l, sin 2 9), and sab 
is the two-dimensional Levi-Civita tensor. Explicitly 
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2 sin 9 cos 9 sin 20, 
2 sin 2 9 cos 20, 



and 



X2 = — 3 sin 2 9 cos 9, 



X] c = cost 

-lc 

v <t> 



Xi c = (2 cos 



sin < 
2 /1 



The tensorial harmonics are defined by Y™ B = K™ B (even 
parity) and X™ B = — X^ A:B ~. (odd parity). Explicitly, 



and 
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= — 4 sin cos sin 0, 

= — sin 2 9 cos 0, 

= — 2 sin 3 #cos#sin0, 
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APPENDIX B: CALCULATION OF THE 
METRIC PERTURBATIONS 

In this Appendix I derive the form of the functions e%, 
62; e3j ^2, and 63 that appear in Sec. Ill B. For this it 
is sufficient to take, say, £12 — £21 and Bx2 = B21 as the 
only nonvanishing components of the tidal fields £ a b and 
B a b- And since the equations for even-parity and odd- 
parity perturbations decouple j^jj, > each case can be 
considered separately. 

Including only even-parity perturbations, Eqs. <|3.5[l - 
(|3.8|) become 

gun = -f(l + f 2 ei£i2sm 2 9sm2(j)), 

2 , - 
Soe = g r e 2^i2 sin cos sin 20, 

2 

Su$ = gf 3 e 2 £i2sin 2 6>cos20, 

gee = f 2 - if 4 e 3 fi 2 (l + cos 2 0) sin 20, 
2 

Se0 = --f 4 e 3 £i2sin6»cos6»cos20, 

= f 2 sin 2 + -f 4 e 3 fi2sin 2 0(l + cos 2 0)sin2^. 

This metric is then substituted into the vacuum Einstein 
field equations. Computing the Einstein tensor is simpli- 
fied by linearizing with respect to £12 and discarding its 
derivatives with respect to u: Since the time scale over 
which £ a \, changes is of order 1Z, the ratio between tem- 
poral and spatial derivatives is of order f/TZ and there- 
fore small in the internal zone; the temporal derivatives 
can be consistently neglected. The field equations pro- 
duce ordinary differential equations to be satisfied by the 
functions e\, and e 3 . Those are easily decoupled, 
and demanding that the functions all approach unity as 
r — > 00 and be well-behaved at r = 2m yields the unique 
solutions 

2m 2 

ei(f) = e 2 (f) = /, e 3 (f) = 1 - — r , 

I 

The first stage of the coordinate transformation is 



as was stated in Eq. 13.9fl . 

Switching now to odd-parity perturbations, Eqs. Q3.5|) - 
(|r?.8|> become 

guu /; 
guf 1; 

2 , - - 

Sue = - ^ b 2 B 12 sm 9 cos 2(j), 

2 

g fi = -f 3 b 2 fii2Sm 2 0cos0sm20, 
2 

gee = r 2 + -f 4 b 3 Bi2COs6cos2<j), 

gg$ = - ^r 4 63^12 sin 0(1 + cos 2 0) sin 20, 
2 

£4,4, = f2 sin2 # ~ o ^&3#i2 sm2 # cos 9 cos 20. 
Following the same procedure, we arrive at 

b 2 (f) = /, & 3 (f) = 1, 

as was stated in Eq. EHUj l. 



APPENDIX C: TRANSFORMATION FROM 
EXTERNAL TO INTERNAL COORDINATES 

Our task in this Appendix is to construct the trans- 
formation from the external coordinates (u,rfl a ) to the 
internal coordinates (u,fCl a ). We shall proceed in three 
stages. The first stage of the transformation, (u,rfl a ) — > 
(u / ,r'n' a ), will be seen to remove unwanted terms of or- 
der m/r in g a p, as listed in Eqs. I|3.37[) - H3.39|l . The 
second stage, (u',r'fl' a ) — ► (u" ,r"fl" a ), will remove all 
terms of order m/lZ in g a '/3'- Finally, the third stage 
(u" ,r"fl" a ) — + (u,fCl a ) will produce the desired inter- 
nal coordinates and return the metric in the form of 
Eqs. 



u = u — 2mlnr, 

- ,a = (i + v)^ 

and it affects the metric at orders m/r and mr/lZ 2 . This transformation redefines the radial coordinate — r — > r' — 
r + m — and incorporates in v! the gravitational time delay contributed by the small mass m. After performing the 
coordinate transformation the metric becomes 

gww = -l-r' 2 £'*+0(r' 3 /^ 3 ) 

+ 7? + h lt + r'(4m£'* - 2a a n' a + + hl^Sl' a ) + 0(mr' 2 /II 3 ), 
= -n> a + lr' 2 (£':+B' a *) + 0(r>yn 3 ) 



19 



tail 
l 0a 



f Ltail 



a "00 

tail 
'OaO 

1 



4m 



(£';+B':)+(5 a b -Q' a Q' b )a b 



I Ltail i Ltailo'fc 



n'h 



a "000 ~r "Oab 



Ltailo'b i O' J,tailo'6 

"Or>b" T I'n'tnOh 1 ' 



a'TOb 1 



+ 0(mr' 2 /^ 3 ), 



^a^b^oo' 1 



/ l tail i Ltail 
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+ W + n' b B f ; ) + n;n' 6 M + hlf c n' c ) + n' a (h^ + h$l^ 



0(mr' 2 /TZ 3 



This metric matches g(internal zone) at orders 1, r' 2 /1Z 2 , and m/r', but there is still a mismatch at orders m/TZ and 
mr'/K 2 . 

The second stage of the coordinate transformation is 



hif(u') du< - \ r ' Ihlfiu') + 2hlf(u')n' a + hf b a (u')n' a n /b 



a b 



(u')x' 1 



and it affects the metric at orders m/TZ and mr/TZ 2 . After performing this transformation the metric becomes 
gu»u» = -l-r" 2 £"* +0(r" 3 /TZ 3 ) 



2m 
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tail 



2 ( a a — -hoQa T lL Oa0 



tail i Q ffa 
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3 V 
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"/^tail i LtailrV/c 
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/ Ltail i l tail(->//c 
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+ 0(mr" 2 /ft 3 ). 



To arrive at these expressions we had to involve the relations 



Ltail _ ttail 

du" 00 ~ 000 ' 



du" 



du" 



^tail 



^tail 
l ab 



■tail 



4m£ ab + hftft, 



which are obtained by covariant differentiation of Eq. Q3.280 in the direction of u a . The metric now matches 
g(internal zone) at orders 1, r" 2 /1Z 2 , m/r", and m/TZ, but there is still a mismatch at order mr" /TZ 2 . 
The third and final stage of the coordinate transformation is 
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Ltail i ( L,tail 



(AS 
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H~ i^Qabc Qbca Qcab 

where 

Qahc = 2^o6c + y (zacdB d b + E bcd B d ^J . 

This produces the metric of Eqs. p.4U|) - l|3.42|) . 
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